Inverse targeting - an effective immunization strategy 
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We propose a new method to immunize populations or computer networks against epidemics 
which is more efficient than any method considered before. The novelty of our method resides in 
the way of determining the immunization targets. First we identify those individuals or computers 
that contribute the least to the disease spreading measured through their contribution to the size 
of the largest connected cluster in the social or a computer network. The immunization process 
follows the list of identified individuals or computers in inverse order, immunizing first those which 
are most relevant for the epidemic spreading. We have applied our immunization strategy to several 
model networks and two real networks, the Internet and the collaboration network of high energy 
physicists. We find that our new immunization strategy is in the case of model networks up to 14%, 
and for real networks up to 33% more efficient than immunizing dynamically the most connected 
nodes in a network. Our strategy is also numerically efficient and can therefore be applied to large 
systems. 

PACS numbers: 64.60.ah, 64.60.al, 89.75.Fb 



I. INTRODUCTION 

The threat of global spreading of epidemics, like the 
pandemic flu from 2009, or spreading of computer viruses 
which are endangering the functioning of Internet depen- 
dent facilities, are responsible for the enormous increase 
in public interest on immunization during the last years. 
Much progress has been achieved in understanding epi- 
demic spreading PI , and various models have been de- 
veloped suggesting possible vifays of efficient immuniza- 
tion [31-HO]. However, the search for even more effec- 
tive immunization strategies must be pursued since any 
improvement of immunization efhciency can save human 
lives and resources. In this paper we introduce a novel 
immunization strategy based on inverse targeting, which 
proves to be effective and numerically more efficient than 
any proposed previous one. 

Epidemics can spread in human population through 
networks of social contacts, and viruses can propagate 
on computer networks. We will implement an immuniza- 
tion process on networks by immunizing the nodes which 
represent either people that are to be vaccinated, or com- 
puters that should be equipped with specially developed 
antivirus programs. 

We suppose that the dynamics of both, the epidemic 
spreading and the immunization process are much faster 
than the growth and change of the network itself. The 
question we want to answer is, given the topology of the 
network through which the epidemics can spread, what 
is the best possible way to immunize its nodes. The im- 
munization process should require the smallest possible 
number of immunization doses. At the same time, the 
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size of the largest connected cluster of non-immunized 
nodes should stay as small as possible throughout the 
immunization process. 

An immunization strategy determines the sequence in 
which nodes of a network should be immunized to pre- 
vent most efficiently the spreading of epidemics. The 
simplest immunization strategy is the one in which the 
nodes are chosen at random |8j. This strategy has how- 
ever proven to be inefficient, requiring a large fraction 
of nodes (in many cases, such as Internet, nearly 100%) 
to be immunized in order to stop epidemics spreading 
[m 112) . Targeted immunization strategies turn out to 
be much more efficient d3j. Here one immunizes first 
those individuals that are most important for the disease 
spreading. In terms of network properties this impor- 
tance can be expressed through the number of links of a 
node, or through its betweenness centrality [131 [TS]. Al- 
though it has been shown that these two properties are 
correlated [7], their role in the immunization process is 
typically quite different. Most studied targeted strategies 
are based on these two node traits. 

In the high degree based (HD)[ni[in], and the high be- 
tweenness based (HB) [T targeting the sequence in which 
nodes will be immunized is based on their degree or be- 
tweenness in the initial network. They can be improved 
by adaptive strategies |7], high degree adaptive (HDA), 
and high betweenness adaptive (HBA) strategy. Here the 
degree or betweenness of nodes is recalculated for the re- 
maining network of non-immunized nodes each time a 
node has been immunized. The node with highest value 
of this recalculated degree (for HDA) or betweenness (for 
HBA) is immunized next. These dynamic strategies have 
proven to be very effective [7]. 

Other strategies, like immunizing neighbors of a ran- 
domly chosen node according to some rule [31 H] have the 
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FIG. 1. (Colors online) Illustration of the inverse targeting on a very small network. In each step of the reconstruction process 
(from left to right) a black node is chosen which when colored green together with its links, gives the smallest contribution to 
the size of the largest green cluster. The inverse of the sequence of nodes obtained from the reconstruction process defines the 
immunization strategy. In each step of the immunization process (from right to left) a node will be immunized following this 
inverse sequence. 



advantage of requiring only local information, but are 
also less effective. A further approach, recently intro- 
duced in Ref. [5] , is the improvement of existing strate- 
gies through optimization procedures. This approach 
leads to very effective node and link immunization strate- 
gies. Another strategy based on equal graph partitioning 
[B^ is especially effective when the number of available 
immunization doses is fixed. For any given number of 
nodes this strategy finds an effective way to partition the 
network in clusters such that the size of the largest clus- 
ter is the smallest possible. In Ref. [B] the authors show 
that for a given number of doses this strategy is more 
effective than targeting HD, HB or HDA strategy. How- 
ever, if the supply of immunization doses is not known 
in advance or increases with time, one must take into 
account that when new doses become available the net- 
work of non-immunized nodes can be quite different from 
the initial network. The nodes which can now be immu- 
nized are then different from those that would give an 
ideal graph partitioning if we had this larger quantity of 
immunization doses from the beginning of the process. 

In such cases dynamic HDA remains the best known 
immunization strategy because the HBA, which needs 
more global information, turn out to be numerically ex- 
tremely demanding. As a consequence, although the 
HBA may be theoretically the more efficient strategy [7] , 
the complexity of its algorithm makes it uninteresting for 
practical purposes and limits its use to very small sys- 
tems. In this paper we will introduce a strategy which 
is more efficient than HDA and in many cases even bet- 
ter than HBA, being at the same time numerically much 
faster than HBA allowing to apply it to large, real sys- 
tems. 



II. INVERSE TARGETING 

In any successful immunization strategy, nodes which 
are not relevant for keeping the largest cluster of non- 
immunized nodes together, should only be immunized 
during later stages. In our method we start by recog- 



nizing such nodes first. This way we reveal the sequence 
of nodes that are going to be immunized in reversed or- 
der. To illustrate our method, let the nodes and the 
links of the network be black. Every node has an index 
(ID) (Figjl]). At each step we check for each black node 
how would turning it and its links green increase the size 
of the largest connected green cluster and then choose 
the one with the smallest contribution to its size. When 
more nodes give the same contribution, the one with the 
smallest number of links connecting it to black nodes is 
chosen. If there are several nodes which fulfill both cri- 
teria, one of them is chosen at random. After a node is 
chosen we put its ID in a list and continue the process 
until all nodes and links are colored green. The list of 
nodes' IDs describing the reconstruction process gives us 
the inverse sequence of nodes that are to be immunized. 

The complexity of the inverse targeting algorithm is at 
most 0{MN), or 0{N'^) for sparse networks, where M 
is the number of links and N the number of nodes. This 
is much faster then the HBA strategy. The best known 
algorithm for determining betweenness of nodes goes as 
0{MN) [ini |T7] and since it has to be applied at each 
step of the immunization process, the complexity of the 
algorithm for the HBA strategy is 0{MN^), or 0{N^) 
for sparse networks. 



III. RESULTS 

We compare the effectiveness of our inverse target- 
ing immunization strategy to HDA and HBA strategies. 
We start by comparing the fractions qc of immunized 
nodes for which the network of non-immunized individ- 
uals breaks apart, which is the measure usually used in 
studies of epidemics. In addition, we use network's ro- 
bustness to immunization as the measure of how effective 
the immunization procedure is. A good immunization 
strategy should make such a robustness as small as pos- 
sible. The robustness sums up the sizes of the largest 
connected clusters S{Q) of the networks of non immu- 
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FIG. 2. The fraction s{q) of sites belonging to the largest connected cluster versus the fraction q = Q/N of immunized nodes 
using the new immunization strategy (green dashed lines) , the HDA attacks (red doted lines) and HDA (black full lines) for (a) 
random regular networks with TV = 8000, M = 16000, (b) Erdos-Renyi networks with TV = 8000, M = 16000 and (c) scale-free 
networks with TV = 8000 and 7 = 2.5 . 
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FIG. 3. The fraction s{q) of sites belonging to the largest 
connected cluster versus the fraction q = Q/N of immunized 
nodes using the new immunization strategy (full lines) and 
HDA attacks (dashed red lines) for (a) the AS Internet with 
TV = 18124 and M = 37357 and (b) the HEP network with 
TV = 27240 and M = 341923 . 
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where TV is the size of the network. This measure cap- 
tures the network response to immunization throughout 
the immunization process, and not only at the percola- 
tion point at which the network of non immunized nodes 
becomes disconnected [18j [T^. The process of immu- 
nization usually takes some time. The initially slow pro- 
duction of needed vaccine soon after the disease starts 
spreading, slows down the immunization process in these 



initial phases. Therefore it is important that the immu- 
nization process not only completely breaks up the net- 
work of contacts among susceptible, non-immunized in- 
dividuals, but also that the size of the largest connected 
cluster, that is, the number of individuals being at high- 
est risk of getting infected, is as small as possible during 
the whole immunization process. This property is ex- 
pressed by a small robustness R. 

We have applied our strategy together with HDA and 
HBA strategies to three different model networks, so 
called random regular graphs ^\ (in which nodes are 
connected randomly keeping the number of links per node 
constant), Erdos-Renyi networks ^21 and scale free net- 
works created with the configuration model |22j . In Fig- 
ure |2] we show how the fraction of nodes in the largest 
connected non-immunized cluster s(g) changes with the 
fraction q of immunized nodes during different immuniza- 
tion processes on these three model networks, q^. is the 
value for which the largest connected cluster is of order 
1/TV, and the robustness R is the area under the curves 
in Figure [2j 

We can see that for all three types of model networks 
our immunization strategy has a smaller qc than HDA 
and has approximately the same value as for the HBA 
immunization. At the same time the robustness of net- 
works in the case of our new strategy is lower than in 
other strategies making it the most efficient of the three 
strategies. 

We have also applied our strategy to two real networks, 
the Internet at the level of autonomous system [53], as 
example for a network through which computer viruses 
can spread, and the collaboration network of high en- 
ergy physicist (HEP) as an example of social networks 
[53]. Here we could compare our strategy only to the 
HDA immunization strategy since HBA would take too 
much computer time. The results are presented in Figure 
[3] where we can see that in both cases the new immu- 
nization strategy is according to both criteria, qc and R, 
much more efficient than the HDA strategy. 

Summing up the results presented in Figures [2] and 
[3] we have found that our immunization strategy out- 
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FIG. 4. The robustness versus the system size TV for random 
regular networks with (k) = 4, Erdos-Renyi networks with 
(k) = 4, and scale-free networks with 7 — 2.5 and femin = 2. 



performs the other two. Comparing to the HDA (HBA) 
strategy, in the case of the scale free network, Qc is 17(0)% 
and R 11(4)% smaller, for the Erdos-Renyi network, Qc is 
24(0)% and R 11(7)% smaller, and for the random regu- 
lar network qc is 16(2)% and R 14(8)% smaller. For real 
networks, for HEP network, is 18% and R 16% smaller 
and for the AS Internet, qc is 37% and R 33% smaller. 

For model networks we have studied the effect of the 
system size on the efficiency of our immunization strat- 
egy and have compared it with the other two strategies. 
In Figure |4] we can see that our immunization strategy is 
for all network sizes more effective than HDA (gives the 
smallest R value). For Erdos-Renyi and scale free net- 
works it is always more efficient than the HBA strategy 
and only for small regular networks with less than 1000 
nodes HBA is more efficient. 

The way the effectiveness of immunization strategies 
changes with the density of links in the model networks 
is shown in Figure [5j We can see that the new immuniza- 
tion strategy is the most efficient for all average degrees 
(k) studied. In all cases studied here (for different (fc) and 
different N) the difference in the effectiveness of different 
strategies is smallest for scale free networks. 



be the most effective immunization strategy [?]■ Due 
to its numerical inefficiency, the HBA strategy can not 
be used for larger systems and is therefore not applica- 
ble to real world problems. Our strategy is in contrast 
also numerically efficient and can therefore become the 
strategy of choice for real applications. The algorithm 
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FIG. 5. The robustness R for all three strategies versus aver- 
age degree (k) for (a) random regular, (b) Erdos- Renyi, and 
(c) scale-free networks with = 8000. 



we introduced could find other applications besides im- 
munization. It could for example be used to determine 
attack strategies for the efficient destruction of criminal 
networks. 



IV. CONCLUSIONS 

We have introduced a novel immunization strategy 
based on inverse targeting. To estimate the efficiency 
of our strategy we have not only used the percolation 
transition qc but also an additional measure taking into 
account the history of the immunization process. We 
find that the new strategy is much more efficient than 
the high degree based adaptive targeted strategy. It is 
also at least as efficient as the high betweenness based 
targeting strategy which was for a long time claimed to 
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